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ABSTRACT 


The study represents a continuation of the development of 
a model output statistics scheme to specify marine fog over 
the open ocean and in coastal waters. Thirty-seven direct 
and derived Fleet Numerical Weather Central model output para- 
meters, monthly climatological fog frequencies, combinations 
of the aforementioned parameters (i.e., interactive parameters) 
and a persistence parameter are used as predictors in a step- 
wise multiple linear regression approach to estimate a predic-— 
tand defined as marine fog probability. The predictand is 
categorized in two ways, in one case (FOGCAT I) as smoothed 
probabilities from 0 to 100% as a function of present weather, 
past weather, visibility and low cloud type; and, in another 
case (FOGCAT II) as a limited number of discrete probabilities 
(to include O and 100%) derived from present weather, past 
weather and visibility only. This study derives diagnostic 
regression equations only using as a dependent data sample 
over 24,000 North Pacific Ocean (30-60N) surface synoptic ship 
observations at 0000 GMT for June through August 1976 and 1977. 
The predictor parameters contributing most significantly to 
the variance are sensible and evaporative heat fluxes, monthly 
climatological fog frequencies, and meridional wind sveed. 
Threat, Heidke skill, and Panofsky-Brier probability scoring 
methods are applied to a selection of the derived equations. 
Predictand variance explained reaches .670, threat/skill 
scores reach 0.42 and probability scores are as low as 0.28 
using the FOGCAT I categorization scheme for the predictand. 
Equations for June and July appear more stable than those for 


August. 
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I. INTRODUCTION 


Marine fog presents a threat to safe nautical activities 
as well as to low-level aviation over the oceans. Histori- 
cal monetary and human losses associated with United States 
Navy operations, attributable solely or mostly to fog, have 
been documented by Wheeler and Leipper (1974). These types 
of losses should be significantly reduced as the methods for 
analysis and forecasting of marine fog become more accurate. 
The environmental science community would additionally bene- 
fit from the increased accuracy of fog analyses/forecasts by 
a commensurate improvement in specifying fog-associated para- 
meters in various analysis models, especially boundary layer 
models. 

Research in the past several years conducted by the De- 
partments of Meteorology and Oceanography, Naval Postgraduate 
School (NPS), Monterey, California, has been directed at more 
adequate diagnoses of marine fog and a better representation 
of its climatology, as based on ocean-station-vessel and 
ship-of-opportunity observations (Renard, Englebretson and 
Daughenbaugh, 1975; Renard, 1976). Additionally, more recent 
research has attempted to use geostationary weather satellite 
data to diagnose marine fog areas (Ihli and Renard, 1977: 
McNab, 1979). The totality and continuity of such satellite 


Observations makes this approach potentially very useful. 


lui) 





However, much further testing is needed before operational 
application is feasible. 

In a more Statistical vein, an earlier study by Nelson 
(1972) developed regression equations uSing the "perfect 
prognostic approach" to forecast visibility (and hence fog) 
at ocean station vessels in the North Atlantic Ocean. This 
approach is not unlike the Model Output Statistics (MOS) 
methods used currently by the National Weather Service (NWS) 
to forecast weather parameters of interest over the conti- 
nental U.S., Hawaii and Alaska (Glahn and Lowry, 1972), ex- 
cept the "imperfect prognostic approach" is the present mode. 
Nelson found quite high correlations of visibility with wind 
speed, relative moisture content and various evaporative 
parameters, uSing observed data. However, when he attempted 
to forecast fog using output parameters from the United 
States Navy Fleet Numerical Weather Central's (FNWC) numeri- 
cal prediction models, his regression equations showed little 
Skill. Earlier work, conducted by Schramm (1966), found high 
correlations between observed values of humidity, air tempera- 
ture, wind speed and visibility at sea; however, the regres- 
Sion equation developed from these parameters was of little value 
at times of low visibility. 

Additional research into the feasibility of developing 
regression equations to forecast marine fog, based on FNWC's 
model output parameters, has been conducted by Van Orman and 
Renard (1977) and Quinn (1978). Both studies developed 


Schemes which statistically processed up to 37 direct or 
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derived FNWC model output analysis—time- parameters in con- 
junction with surface synoptic ship observations to 

generate multiple linear regression equations yielding as a 
predictand the probability of marine fog occurrence at 
analysis time over the North Pacific Ocean. The Van Orman 
/Renard results, based solely on July 1976 data, showed 
considerable skill over FNWC's advective fog model, FTER, 

as applied to the analysis of fog (U.S. Naval Weather Service 
Command, 1975) and climatological fog frequencies developed 
at the NPS (Willms, 1975). Quinn expanded the data base to 
also include June and August 1976 and developed regression 
equations both for the unique probabilistic marine fog pre- 
dictand introduced by Van Orman and Renard as well as a 
modification of that predictand. Quinn further introduced 
interactive predictor parameters into the regression scheme 
by combining the climatological fog frequency parameter with 
the two most Significant model output parameters to derive a 
new set of diagnostic regression equations. A clear improve- 
ment over climatology and FNWC's FTER was demonstrated using 
these modified regression equations on dependent data. The 
equations based on the June 1976 data also showed improvement 
over FTER and climatology when used on an independent data 
Set from July 1976. These prior studies give clear indication 
of the relative accuracy of the MOS approach in diagnosing 
marine fog, with a successful extension to forecast modes yet 


to be demonstrated. 


Some parameters are taken from analysis models, other are 
diagnostic parameters from prognostic models. 


Ie 








Currently, only FNWC is producing fog analysis/forecasts 
On a large scale over the oceanic areas through their sta- 
Pst recal prebDabllistic product, FTER. A climatological 
parameter is not used in their product due to the unavail- 
ape y OL an accurate climatology at the time of its 
development. A comprehensive uniform climatology of fog 
occurrence over all of the Northern Hemisphere is now avail- 
able from the National Climatic Center (Guttman, 1978) and 
has been incorporated into this study. This inclusion re- 
presents a refinement of the fog parameter from the two 


prior studies. 
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II. OBJECTIVES AND APPROACH 


The primary objective of this study was to continue the 
development of a multiple linear regression approach begun 
by Van Orman and Renard (1977) and pursued by Quinn (1978) 
to specify effectively the distribution and likelihood of 
Marine fog at analysis time over the middle latitude 0-— 
60°N) region of the North Pacific Ocean during the summer 
season. Part of the primary objective was to evaluate the 
equations’ skill in comparison to the climatological fog fre- 
quencies derived by Guttman (1978) and FNWC's existing opera- 
tional scheme for forecasting advection fog probability, 
FTER. A second objective was to determine whether the pre- 
dictand categorization scheme FOGCAT I developed by Van Orman 
and Renard or the FOGCAT II scheme developed by Quinn is 
Sigetior In declining the occurrence of fog. Additionally, 
an investigation of the usefulness of persistence as a 
diagnostic prediction parameter was initiated. The approach 
follows that of Quinn's (1978) experiments with the addition 
of June, July and August 1977 data. 

Sets of regression equations were developed for both fog 
mrvevorization schemes for (a) each of the six months June, July 
end August 1976 and 1977; (b) the combined data for both 
Junes, both Julys and both Augusts, (c) the combined data 
for June, July and August 1976, (d) the combined data for 


June, July, and August 1977, and (e) all of the data for the 
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Six-month period. Those equations which were based on data 
from both June's, both July's, both August's, and all six 
months were scored using three methods, skill, threat and 


provabLllity. 


—ow 
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TIIil. DATA 


A. AREA 


The area of study was confined to the North Pacific Ocean 
to maintain continuity and consistency with the previous 
studies of Van Orman and Renard (1977) and Quinn (1978), 
specifically the area north of 30°N and south of 60°N (Figure 
1). A 27x14 grid was superimposed on this region with grid 
points coincident with their counterparts on the same subarea 


of the standard FNWC 63x63 grid (Figure 2). 


B. TIME PERIOD 

Due to the high frequency of marine fog in the summer 
months over the North Pacific Ocean (Willms, 1975; Guttman, 
1978), the additional summer months (June, July and August) 
of 1977 were added to the available set of the same months 
gm 1976. At OOOO GMT the entire North Pacific Ocean is in 
daylight with local noon occurring at the international date- 
line. Only data from this time were used since the aAecur acy 
©f visibility reports from ships transitting the area should 


be best during daylight hours. 


C. SYNOPTIC WEATHER REPORTS 

The June-August 1976 and 1977 synoptic weather reports 
used in this study were provided by the Naval Weather Service 
Detachment, physically located with the National Oceanic and 


Atmospheric Administration's National Climatic Center, 


IEG 





Asheville, North Carolina. These data had been edited to 
eliminate duplicate reports. The total number of observa- 
tions received were 4277 for June 1976, 4391 for July 1976, 
4134 for August 1976, 4232 for June 1977, 4198 for July 

1977 and 3156 for August 1977. All observations were ob- 
tained from transitting ships and Ocean Weather Station 

Papa except data obtained from 11 land stations located in 
the study area (Appendix A). As all of these land stations 
are located on islands or immediately adjacent to the coast 
and have relatively low elevations, the effect of topography 


is thought to be minimum. 


D. MODEL OUTPUT PARAMETERS 

FNWC provided the 22 diagnostic model Output parameters 
(MOP's) for the time period and area of interest as Output 
from several of their numerical analysis and prediction 
models: the Mass Structure Analysis Model, the Primitive 
Equation Prediction Model, the Marine Wind Model, and the 
Spectral Ocean Wave Model (U.S. Naval Weather Service Command, 
1975). Additionally, a set of 15 parameters were derived 


from the primary set and eight interactive parameters were 


developed. A complete listing and description of all these 


Darameters 1S located in Appendix B. 


E. FOG FREQUENCY CLIMATOLOGY 
Previous research used the Naval Postgraduate School's 
North Pacific Ocean marine fog climatology (Willms, 1975). 


AS an updated hemispheric fog frequency climatology (Gittman, 
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1978) is now available, it was used exclusively throughout 
this research. Specific values for fog frequency were inter- 
polated to an 89x45 grid in the area of study, which is 
approximately four times as dense as the FNWC grid over the 
region. This finer resolution allowed a more accurate speci- 
meecatlon of the fog frequency climatology. The fog frequency 


climatologies for each month are presented in Figures 3-5. 


IS, 





IV. PROCEDURE 


A. PREDICTAND CATEGORIZATION METHODS 

The two methods of predictand (fog) categorization used 
throughout this study are based on four elements of the 
Synoptic reports, namely: present weather (ww), past weather 
Sipeeevesibility (VIS), and low cloud type (CL). Observation 
elements other than present weather are used due to incon- 
Sistencies in the reported observations, which result in 
part trom current rules for coding observed data (U.S. De- 
Pomements Of Commerce, Defense and Transportation, 1969) 
(See Appendix C). For example, fog is not coded as present 
weather whenever any form of precipitation is occurring 
Simultaneously. Both methods of fog categorization assign 
a probability of fog occurrence to each synoptic report, 
differing mainly by the number and combination of observa- 
tion elements used to assign the probability value. 

The Fog Categorization I (FOGCAT I) scheme was developed 
by Van Orman and Renard (1977) and used all of the four 
elements listed above. This scheme assigns one of five 
major fog categories and one of three subcategories to each 
Beymoptilc report as a function of ww, W, and CL. Then a fog 
probability, O to 100% is assigned to each observation de- 
pending upon the major category and subcategory assigned. 


Appendix D outlines FOGCAT I. 


20 





The Fog Categorization II (FOGCAT II) scheme developed 
by Quinn (1978) uses only present weather, past weather and 
Moi iityaromssien one or Seven discrete fog probabilities 
(0%, 10%, 15%, 35%, 60%, 85% and 100%) to each synoptic 
Observation. This scheme is more diréct in its method of 
categorization and assigns fog probabilities based on a 
filtering technique. First present weather is evaluated. 
If present weather definitely indicates either the presence 
or absence of fog, then a probability of 100% or 0% is 
asSigned, respectively. In either of these cases, both past 
weather and visibility are ignored. For all other present 
weather codes, past weather and visibility are considered 
in the assignment of a fog probability. Past weather is 
evaluated first. If the possibility of fog is eliminated by 
the past weather code, a score of 0% is assigned and 
visibility is not considered. If the past weather code in- 
dicates the possibility of fog, then the visibility is taken 
into account in assigning one of the choices for an inter- 
Meclate Loe probability (i.e. 10, 15, 35, 60, 80). The 
asSignment of these intermediate values 1S somewhat arbitrary. 
Details of FOGCAT [II are presented in Appendix E. 

Neither the FOGCAT I nor FOGCAT II schemes represent an 
ideal categorization scheme for fog probability. FOGCAT I, 
on the one hand, could assign an observation a fog probability 
aS low as 62.1% even though the present weather code definitely 
indicates the presence of fog, or it could assign a proba- 


bility as high as 24.1% when fog appears very improbable from 


Pah 
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the synoptic report. FOGCAT II was developed to achieve a 
Bemrer fit Of probability to occurrence of event. Except 
for the "100" and "0" values, the FOGCAT II probabilities, 
like those of FOGCAT I, are somewhat arbitrary. In any case 
the schemes were accepted as developed by the previous 
authors, recognizing that some adjustments may improve the 
physical representation of the categorization along with its 


verification. 


B. PREDICTOR INTERPOLATION 

The interpolation method used to determine values of the 
predictors at the observation points iS a natural bicubic 
spline curvilinear interpolation scheme. The locally de- 
veloped program, SPLIN, is available at the NPS W. R. Church 


Computer Center and was used throughout the study. 


C. REGRESSION SCHEME 

A stepwise multiple linear regression program, called 
BMDP2R (University of California, 1975) was chosen as the 
means for deriving the marine fog diagnostic equations. The 
polynomial regression scheme also availabe in the BMD com- 
putCY program series was not used due to its inability to 
handle data sets of over 1000 cases. As each month has 
approximately 4000 synoptic observations and combined data 
Sets had at least 7000 cases, the usefulness of the polynomial 
regression scheme applied to smaller segments of the data was 
unlikely to give accurate results. A non-linear, non-poly- 


nomial scheme alSo available in the BMD series uses excessive 
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computer memory for a 45-variable regression, and, therefore, 
was likewise not used. 

BMDP2R computes a Sequence of multiple linear equations 
in a stepwise manner. At each step one variable is either 
added or removed from the previous step'’s equation as dependent 
on the F-to-enter and F-to-remove criteria. In the forward 
selection procedure the dependent predictor inserted in the 
equation is the one with the highest coefficient of partial 
fectermination (AR“). At each step the regression procedure 
meevaluates the variables already in the equation, and may 
find that a variable important at an earlier stage may be 
less important at a later stage due to the intercorrelation 
of variables. In this procedure, before a variable is added, 
the variable already in the equation with the partial F 
value is dropped if this latter value is less than the maxi- 
mum F-to-remove value (Wesolowsky, 1976). Selection proceeds 
until no candidate variables qualify or until all independent 
Variables have been used. 

The BMDP2R program permits the specification of the 
tolerance level. If the tolerance is near zero, the variable 
being considered is close to being a linear combination of 
the variables already in the equation. Variables with low 
Tolerances could cause computational difficulties in the 
method of calculation (Wesolowsky, 1976). Since tolerance is 
defined as one minus the coefficient of determination, the 
value of .0Ol1 was used as the minimum acceptable tolerance 


menberion (University of California, 1975). 
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The approach taken in this study was to treat the y-inter- 
cept (Y-INTCP) as a variable and to use the BMDP2R program's 


default values of 4.0 and 3.9 as the values for the minimum 


F-to-enter and maximum F-to-remove criteria. 


D. VERIFICATION SCORING 


Two types of verification scores are used in the study to 
test the skill of different regression equations in describing 
the distribution of fog probabilities. The first type, 
exemplified by the Heidke Skill Score (HSS) and the Threat 
Score (TS), is used to test the effectiveness in specifying 
discrete occurrences of fog. The second type (Panofsky-Brier 
Probability Score) tests the accuracy of probabilities in 
estimating the likelihood of fog occurrence (PS) (Panofsky 
and Brier, 1958). 

The formulae used for the three scores (HSS, TS, and PS) 
are given in Appendix F. Currently the Threat Score is 
popularly used by the Techniques Development Laboratory, 
National Weather Service (Bermowitz and Best, 1978; Miller and 
Best, 1978), while the Heidke Skill Score has been traditional- 
ly used by meteorologists. The Panofsky-—Brier Probability 
Score is suitable for evaluating the effectiveness of a 


predictor given in terms of a probability (0-100%). 
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V. RESULTS 


A. REGRESSION EQUATIONS 

The sets of equations generated for the individual months 
of June, July and August 1976 and 1977 as well as combina- 
tions thereof, uSing categorization schemes FOGCAT I and 
FOGCAT II, are listed in Tables I-XII. They are shown in 
stepwise order and include Y-intercept, regression coeffi- 
Clients and the amount of variance explained with the inclu- 
Sion or removal of each variable (R? and AR@). Only six 
steps are given in each table for purposes of brevity. In 
all cases no two successive steps added a cumulative AR? Of 
greater than .005 after the sixth step until a variable was 
removed. Recognizing that variables which explain less than 
1% additional variance are normally not included in opera- 
tional regression equations, six steps were presented so 
that all variables having a possible effect on further devel- 
Opment could be identified. 

Since the climatological synoptic regime varies from 
month-to-month and from year-to-year the regression equa- 
tions vary accordingly. Several Similarities are evident, 
however. First, in all cases for both FOGCAT I and FOGCAT II 
the first model output parameter to be entered is one of the 
heat flux terms, either EHF or SEHF. Second, either CLIMO 
Or an interactive parameter containing CLIMO enters all the 


equations by step 3 with the exception of the equations 
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derived from August 1977 data and the equation for all 
months for FOGCAT I. Finally, the meridional component of 
the wind, v, or an interactive parameter containing v enters 
all of the equations no later than step 4 with the exception 
of the equations derived using FOGCAT I on the combined 

Moy O7G/7?% data or FOGCAT I and II, on the August 1976 

and FOGCAT I on the August 1977 data. Additional similari- 
ties are evident at later steps but are not considered 
Sreniticant tO this discussion. 

In all cases, the first three non-constant variables 
entered into the regression equations accounts for between 
71.9% and 92.5% of the summation of coefficients of partial 
determination CAR“) determined by the entire set of model 
Sacpuc predictors. 

As the magnitudes of the coefficients of determination 
Cie) vary greatly between those equations derived using the 
two different categorization schemes, a direct comparison 
of the effectiveness of the equations based on FOGCAT I 
versus those based on FOGCAT II is difficult. Although Table 


XIIE indicates that the ZAR“ 


for FOGCAT II is generally 
greater than this quantity for FOGCAT I, comparisons are more 
properly made using standard verification scoring methods, 
mein section B following. 

The eight interactive parameters used in the regression 
analysis were formed by either multiplying the parameters 
found to be most significant by Quinn (1978) by climatology 


Or by multiplying various of these most significant parameters 


by other significant parameters. The inclusion of these 
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interactive parameters allows climatology to enter the scheme 
at an earlier step. Originally, climatology was too highly 
correlated with the heat flux terms to enter the equations 
alone, and thus it was rejected by the tolerance criterion. 
Having identified the most probable meteorological indi- 
cators of fog occurrence by regression analysis, the feasi- 
bility of adding a persistence parameter was investigated. 
Due to lack of time, this investigation was limited to the 
iy Logo and 1977 data sets. The persistence parameter 
(PERS) was based only on those observations which had at 
least one synoptic report within a one degree latitude/longi- 
tude box of the reported position on the previous day. This 
requirement reduced the data set for July 1976 from 4391 to 
36a reports and for July 1977 from 4198 to 1591 reports. 
Regression analysis, with PERS included, was then per- 
formed on the July data sets, both individually and in unison, 
uSing both categorization schemes. The sets of equations for 
these data sets are listed in Tables XIV-XVI. Persistence 
displayed a high correlation with EHF, SEHF and CLIMO and 
was rejected for inclusion in the equations by the tolerance 
criterion on this basis until the seventh step at the earliest. 
The Similarities among equations which were noted earlier are 
decreased when PERS is included with only the heat flux term 
metalning its former importance in all cases. This instability 
introduced into the equations may be only the result of the 
reduced data set; the investigation of persistence as a fog 


indicator should not be abandoned until it is tested on a 


larger scale, both alone and in interactive combinations. 
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Fields of selected model output parameters, FTER and re- 
sression-computed fog probabilities for the arbitrarily 
selected observation time, 0000 GMT 3 August 1976, are shown 
in Figures 6-12. Considerable frontal activity in the mid- 
latitude central and western North Pacific Ocean is indicated 
by the sea-level pressure analysis (Figure 6). The two most 
important fog-related parameters, EHF and SEHF, are shown in 
Figures 7 and 8. Maximum positive values at the subtropical 
latitudes give way to a broad zonal band of negative and low 
positive values centered between 45-55N. Another important 
regression parameter, the algebraic value of meridional wind 
Speed, v, is displayed in Figure 9. This field closely re- 
lates to the geostrophic/gradient wind implied by the sea- 
level pressure analysis. Figures 10 and 11 show the regres- 
Sion-computed fog probabilities based on three variables for 
Fog Categorization Schemes I and II. The configuration of 
these fields resembles that of SEHF, the FNWC variable ex- 
plaining most of the variance in fog probability. Probabili- 
ties on the FOGCAT I version (Figure 10) appear to be about 
15% greater than those of FOGCAT II (Figure 11). The differ- 
ence 1S particularly noticeable in the subtropics where fog 
mecelihood 1S climatologically small. The FTER probability 
field (Figure 12) shows a smaller range of values than the 
NPS probabilities, and is a closer fit to the FOGCAT II 


regression probabilities. 


B. VERIFICATION SCORING 
For purposes of verifying fog occurrence in an observation, 


major categories S and F were used in the FOGCAT I scheme, 
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Bnd an assigned probability of fog occurrence of 100% was used 
in the FOGCAT II scheme (Appendices D and E). An attempt was 
made to use 85% fog likelihood as also verifying fog occur- 
rence in FOGCAT II, but this criterion yielded significant- 

ly poorer results than 100% used alone, as had been found 

by Quinn (1978). This may indicate that some fine tuning of 
the assigned fog probabilities in FOGCAT II would yield 

Better results. 

Optimum threshold probabilities were also computed for 
those equations based on the combined June data, combined 
July data, combined August data, and combined data for June, 
July and August 1976 and L977. The threshold probability is 
that equation-specified probability which best separates the 
fog/no fog occurrences, hence yielding the best skill and 
threat scores (Bermowitz and Best, 1978). The threshold 
Heidke Skill Scores and Threat Scores presented in Tables 
AVIIT-XXTI are those which maximize the indicated scores. The 
Panofsky-Brier Probability Score presented in these tables 

‘directly evaluates the accuracy of regression probabilities 
| in estimating the likelihood of the event. 

The regression equations as derived by the BMDP2R program 
were not used in their entirety in the verification phase of 

| this Study nor were all the equations verified. Only the 
Seat ions in the form containing the Y-intercept and the 
first three non-constant variables were used. It is possible 
to "overfit" the regression equation to the dependent data. 


Therefore only variables which bear a physical relationship 
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to the predictand should be included (Panofsky and Brier, 
1958). In all the equations verified, the first three vari- 
able predictors included had an apparently strong physical 
relationship with the occurrence of fog, while several equa- 
tions introduced terms with only a weak physical relation- 
ship to fog likelihood at step 5. 

The verification scores for the five data sets for which 
regression equations were evaluated and compared to ove veri-~ 
fication scores for FTER and climatology are presented in 
Tables XVII-~XXI. Note that the scores for the regression 
equations as well as for FTER and climatology vary from the 
FOGCAT I to the FOGCAT II schemes. Since verification in 
each scheme is against the actual occurrence of fog (Category 
S or F in FOGCAT I or 100% in FOGCAT II), a relative compari- 
son of the effectiveness in categorizing fog occurrence 
between FOGCAT I and II for purposes of regression analysis 


is possible. This is not to imply that observations verify 


| as fog/no fog equally with FOGCAT I and FOGCAT II; there are 


| 


Slight differences not considered important in this analysis. 


In four of the five cases verified, the skill and threat 


scores for the regression equations were higher for FOGCAT I 


| than for FOGCAT II. In all verification scores, the FOGCAT 





Il regression equation outperforms both climatology and FTER, 
except for the regression equation based on the two August 


data sets. The FOGCAT I regression equation outperforms 


Climatology and FTER in skill and threat scores for each data 


set except the full June, July and August 1976 and 1977 set 
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where climatology is the top performer. In probability scor- 
ing, FTER outperforms the FOGCAT I regression equation in 
four ok ie five cases and climatology outperforms it in the 
other cases in the P-score. 

The variations of skill and threat scores for climatology 
from FOGCAT I to FOGCAT II may be taken as a measure of the 
differences in the verifying scheme (i.e. strong fogger, 
category S, and foggers, category F, for FOGCAT I and observa- 
tions assigned 100% probability in FOGCAT II). Thus, a mea- 
sure of the relative comparison of FOGCAT I and II regression 
equations might be obtained by noting the changes in the 
relative improvement over climatology. For example with 
reference to Table XVII, FOGCAT I skill/threat scores exceed 
climatologies skill/threat scores by .107/.072, while for 
mOoGCAT if the numbers are only .082/.063, indicating greater 
merit for FOGCAT I. Table XXII tabularizes this information. 
On this basis relative improvement over climatology is best 
} uSing FOGCAT I for equations in Tables XVII, XVIII, XIX 
(except for threat score) and XXI, while FOGCAT II is better 
for equations in Table XX. 

When the persistence parameter is included in the regres- 
Sion scheme for the two Julys, regression skill and threat 
Scores generally show improvement (Tables KVIII-XIX). Con- 
Sidering the reduced data set, at best, it is to be regarded 
aS giving encouragement to the future refinement of a more 


comprehensive persistence parameter. 


eal 





It is to be noted again that all regression equations 
apply to analysis time, hence are diagnostic. Further all 
are relative to dependent data. The data were considered 


too limited to expend even one month as an independent set. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


From the evaluation of these diagnostic Model Output Sta- 
tistics regression equations, it can be seen that the MOS 
moproach to prediction of marine fog holds great potential. 
The developed regression equations show improvement over both 
elimatology and the FNWC advective fog forecasting product, 
FTER. However, it remains difficult to recommend, without 
Sualifications, the use of either FOGCAT I or FOGCAT II 
schemes for the predictand. 

The comments that follow relate to the equation for each 
month (for two-year period); the equation for all six months 
combined are not considered to be operationally useful. 

Even though FOGCAT II equations score better than CLIMO and 
BIER (except for Ausust) uSing all the scoring methods, 
threat, skill and probability, FOGCAT I equations show higher 
threat and skill scores than FOGCAT II and excel CLIMO and 
FTER in two of the scoring methods, threat and skill. It 
appears quite significant that FOGCAT I shows most skill 
relative to climatology and FTER as well and therefore this 
predictand categorization appears to be the best to use. 
Further development of FOGCAT I and II or another approach 
may well yield still further improvement in regression results. 

The following recommendations are offered as a guide to 


future work in this area: 
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li In order to Simplify the analysis procedure and the 
amount of data manipulation, parameters which have little or 
no effect on the outcome of the regression equation should 
be deleted. These include: SOLARAD, DDWW, PPW, PDW, SPW, 
SDW, WCP, THETAX, THETAR, STABX, ASTDX, ASTDR, CAPV, ADTSEA, 
AASTDX, AASTDR, and SSTA. See Appendix B. 

2. The development of a more realistic persistence para- 
meter should be continued and included in the regression 
scheme both alone and in interactive combinations. Asa 
first approximation the fog probability assigned by the 
regression equation (without inclusion of a persistence para- 
meter) to each grid point of the previous day could be inter- 
polated to the actual observation positions of the present 
day. This would give a persistence value to each observation 
whether or not a report was available within a one degree 
latitude/longitude box on the previous day. These first 
approximations could then be weighted for those observations 
which had a report within the specified box on the previous 
day. In this manner the size of the data set would not be 
reduced aS it was in this study. 

3. Additional data (perhaps one more year) should be 
added to further stabilize the regression analysis. 

4. A stratification of the set of regression equations 
by latitude, longitude, or by meteorological phenomena (i.e. 
wind direction, positive/negative thermal advection, or 
Places where climatological fog frequencies are greater) is 


likely to improve accuracy. 
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9. The regression approach should be extended to regions 
other than the North Pacific Ocean and time periods other 
than summer. 

6. As FNWC predictive model output parameters fields 
become avallable for testing they should be incorporated into 
the regression equations for forecast intervals to at least 


48 hours. 
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TABLE XIII. Changes in amount of variance explained 
CAR om different fog categorization 


methods. 















Summation of AR? Due to Inclusion of 
First Three Non—Constant Variables 
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086 


43 














900° 6GV- G09" iG ie ewe SiGe ee OSS" 8PL t- | peaowsey| 9 
cop | 9 e- [ oss” | see'te | test 

aa cae P99 608'0c | pec’ s OLS 9EL'T TEs‘ Tm | G 

“ie Waal | ee | ee ee 

=) a ee ee 

| oor TT pore [vaote | 1 
pee | | |. | bk). ae 

ou ol Xa HIVLAV | UGLSV XGLSVV JHYS | dOLNI-A| das 


(sdois XTS) gogo’ = HV T¥IOL 
"bZO° FO zHV TBUOTITPPe ue poutetdxe pue doys 416 Sut 1% poisezue sous STsJod ‘sdoys 
QT 19a1VFe weardoid yzdaqnd AQ poyeUTWsS4 dSutddey4g ‘“sweyos uot zezta1odeyeO TI LVODOT ad 


sG07 1_vs9"_[eso'ot | 260% [eo e-| teo- | eve’ | coo't-_| 6va've 
200° lives | stort ore't-| zeo- | tte’ | tet | 8v6ve 
es poe | veo | see 908° - | voLse 
Z00° ice ee oe LUV | vbL’ ~ | 096°%e 
100" er EE 
SLO TCC 


ELG LE 0 


si Ny ot sty wy © 





(sdeqys XTS) [60° = HV T¥IOL 
"8TO0° 9 2UV TVeuOTIIppe Ue poute,tTdxo pue dois uz, 943 78 porto US s0uU01STSIeq ‘Ssdo41s 
Gh 19177 wers0id yzdand Aq pozyeUTWsS8 4 sSutdds3g “9woeyos uotzyeztzose ye) I LVYODOA "yy 


-peyuesoid sdets 9g Ajuo AytTaAssrq A0g *popn Tour Iloyouvarzed aoue ys tsaed 
a  . §SNUDTOTTToOU UOTSS ouvue astmde 1s "AIX HIAVIL 








2 | aHas | dOLNI-Al das 


49 





ETE SOS ° 9¢4° {LLb'ETs- | OTb'9 | ESb'T |[abvs're- | Bzo's- | O8P'8T 9 
ol oll ULV ddvVaLs HHS Mdd XQGLS VV HHS COLNI-A duds 


(Ssdoqs XIS) 17i 9] suv o TO tomas 
JO ,HV TBUOCTITppe uv pouTe{Tdxs puv do4s YIEFT 241 1% petojJUS soUSeIsSTSJeg ‘sde4s 
ET fase wesdsoid yzdqnd Aq peyeuTWs94 surddei1g ‘oswayog UOTIeZTIOS91%O ITI LvyoOo0da “ad 


 — ae et ee 
<< a a 


20 


map Ny] Oy Si my © 


| B89" (| LEO’ b L€6°9G 
LLO’ 
UV JOLNI-Al dhs 


(sdeis xts) 960° = 2,4V T¥430L “S00° 
JO ,UV TeuUOTITppe ue pourTetTdxs pue dois YIP 944 1% pete US sdUeysSTSJ9q ‘sdoa4s 
PL tease wessoid yzdaqnd Aq peyeutwste4y sutddsig ‘saweyos UOT ZezZtT£LO99189D [ LVOODOT 'V 





*pojyuoserid sdeis g ATuoO AjTA9Ig IOq ‘“‘pepnToUuT soyjouvred sousystsaed ayy 
UTM ‘S3TO0d9eI TEcT ‘B1BP LZLGT Ane ‘SyUuSeTOTZFIOO UuOTSserdaq ostmdeyisg :AX ATEAVE 














sar sor’ | Sb9o'l | tee zzz] ozo’ fogo'o6z-| L00°t- | zz e- | 86°82 9 
coo he Ed tO” ezerott| oo" [eteten [sseze | + 
“Tae Oa CO 
a 519 | 808° b= | 086 GF 
ee Le 
guy Zt Yiveal XAVES |. oA 4) uaVves pK ana | dO4NI-A‘{ daus 


(sdeis XTS) €bT° = zHV T¥IOL “UOT PeUTWIIy STOO iojyueo jou ptp soueistsdJeq ‘*sdois 


QT rea7e weas0id ycdadNnd Aq po .eUuTWISE4 Sutddojig ‘owoyog uotyeztsodse ye) [I LVODOT ‘d 





Ost | 860. Gus = P60 Tc 
et ae 
a; ee 


ol 


my NY OO, ST uy] © 


XA OWTTO 4H dOLNI-A dubs 


(sdeis xts) ¢60° = ,HV T¥IOL ‘OTTO’ 


JO ,uV [TVUCTITppe ue poutetdxe pue deqis Yyapl ey) 7% petoejue soUueysTSdod ‘sdoqs | 
CT Jo1ye wesis0id yzdand Aq po ,euTWse4 Sutddeig ‘owayog uotyeztaodeyeD IT LVODO V 


: ; d aoueystsiod oy. YATM 
ajzuesoid sdois g Atuo AqjtaAosiq IOg popntTour ra ,owVse | 
: ‘sqzzodea 6g6z ‘B21eD LLGT pue Q1GT Aur fsqueToOTZye00 UOTSSerdoI astmdersg ‘LAX ATEVL 


OO 





TABLE XVIL Skill, threat and P-scores June 1976 and 1977 


data. 
9 
A. FOGCAT I Categorization Scheme. R = .627. 
REGRESSION 
EQUATION FTER CLIMATOLOGY 
t 
Skill Score Ono Reiley 2206 
Threshold 40 14 au 
* 
Threat Score Reo re. 2006 . 266 
Threshold 38 dk 16 
* 
P-Score e297 . 209 ges: 


Boe Probability (%) = 36.725 - 2.095 (EHF) 
- .030 (CLIMO)(SEHF) + .226 (CLIMO) 


2 


B. FOGCAT II Categorization Scheme. R= .367. 
REGRESSION 
EQUATION FTER CLIMATOLOGY 
* 
skill Score 02 AOL 2220 
Threshold 28 23 mage 
* 
Threat Score PES Ee 2200 Pee, 
Threshold 26 at 22 


He 
P-Score .288 339 mo Ole 


Fog Probability (%) = 31.154 - 2.956 (EHF) 
- .032 (CLIMO)(SEHF) + .023 (CLIMO)(V) 


* 
Indicates best score in each category 
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TABLE XVIII. Skill, threat and P-scores July 1976 and 1977 


oleh eek: 
A. FOGCAT I Categorization Scheme. R“ = .670. 
REGRESSION 
EQUATION FTER CLIMATOLOGY 
*k 
Skill Score . 384 P330 5S OFS) 
Threshold 46 35 So 
* 
Threat Score 404 Looe . 360 
Threshold 42 23 28 
*k 
P-Score P3a0 w20 6 ae, 


Mog Probability (%) = 34.8398 - 1.663 (EHF) 
~ .026 (CLIMO)(SEHF) + .292 (CLIMO) 


2 


B. FOGCAT II Categorization Scheme. R = .441. 
REGRESS ION 
EQUATION FTER CLIMATOLOGY 
ok 
Skill Score .378 T3035 , 302 
Threshold 38 35 42 
oe 
Threat Score . 386 » 346 2347 
Threshold 34 23 35 
oe 
P—Score 299 RAPA yB22 


Fog Probability (%) = 37.131 — 4.037 (EHF) 
+ ,O29(CLIMO)(V) + .466 (EAIR) 


* 
Indicates best score in each category 
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TABLE XIX. Skill, threat and P-scores July 1976 and 1977 
data with persistence parameter. 


A. FOGCAT I Categorization Scheme. R~ = .673. 





REGRESSION 
EQUATION FTER CLIMATOLOGY 
* 
Skill Score ~405 334 319 
Threshold 48 28 283 
ok 
Threat Score 423 ree veal TOOL 
Threshold 47 Die) oO 
4 
P—-Score 356 A Se Dae 


Fog Probability (%) = 17.869 - 2.945 (EHF) 
+ .510 (CLIMO) + 1.113 (EX) 


B. FOGCAT II Categorization Scheme. R? = ,458. 


REGRESSION 
EQUATION FTER CLIMATOLOGY 
Pee score 394" 337 320 
Threshold 36 S58) 41 
*« 
Threat Score .414 2305 Aero il 
Threshold 36 25 33 
2 
P-~Score moo ,o10 oo 0, 


Fog Probability (%) = 28.167 -~ 3.956 (EHF) 
+ .642 (TX) ~ 130.026 (STABR) 


2K 
Indicates best score in each category. 
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TABLE XX . Skill, threat and P-scores August 1976 and 
Vo¢d. da ca.. 


A, FOGCAT I Categorization Scheme. R” = .660 


REGRESSION 
EQUATION FTER CLIMATOLOGY 
we 
Skill Score ,40'5 374 . 146 
Threshold 45 28 Z6 
 o 
Threat Score HOLS » 345 219 
Threshold 49 23 14 
ok 
P—Score BAL AS, mL ifs paoZ 


Fog Probability (%) = 38.781 - 1.859 (SEHF) 
+ 1,059 (V) - .708 (CAPU) 


B. FOGCAT II Categorization Scheme. R° = .420. 


REGRESS ION 
EQUATION FTER CLIMATOLOGY 
* 
Skill Score 402 Ou plese 
Threshold 46 23 25 
* 
Threat Score 368 . 341 RPAOKS 
Threshold 40 a4 16 
* 
P—~Score AAS) TE 20 sooo 


Soe Probability (%) = 26.352 - 2.105 (SEHF) 
+ .727 (V) - .067 (SEHF)(V) 


* 
Indicates best score in each category. 
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TABLE XXI . Skill, threat and P-scores June, July and August 
I9%G andeon( data, 


A. FOGCAT I Categorization Scheme. Re = ,646. 





REGRESSION 
EQUATION FTER CLIMATOLOGY 
k 
Skill Score .214 . 196 SSS 
Threshold | A EAS 44 
* 
Threat Score HeLOL AOL .319 
Threshold 16 Al 4] 
mK 
P-Score , 304 azo ,208 
Fog Probability (%) = 31.966 - 1.559 (SEHF) 


Toole Vo) = ae (CC lEEMO ) 


ry 


i FOGCAT Il Categorization Scheme. R = .400. 
REGRESSION 
EQUATION FTER CLIMATOLOGY 
mK 
Skill Score 310 wal ~ 294 
Threshold Se 34 4n 
k 
Threat Score 302 wooo mas OL 
Threshold Sal 23 a0 
* 
P—Score 283 322 -318 


fog Probability (%) = 30.406 - 2.811 (SEHF) 
PeeO som Colony) ts le530 CSHE) 


*K 
Indicates best score in each category. 
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tobe hie Comparisons or ckill and Threat Scores for 
the Regression Equations and Climatology 


FOGCAT I FOGCAT II 
Skill Score |Threat Score |Skill Score ;Threat Score 
Period |Difference Difference Difference Difference 
(Regression |(Regression (Regression | (Regression 
Equation - Equation - Equation - Equation - 
Oematorecy) |climatrology) | Clamatology) |Climatolocy ) 
SO 7 SOW 2 .082 .063 
.081 .044 .O076 039 
.086 062 .O74 “053 
S2on Pikes, .266 . 160 
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Figure 2. Fleet Numerical Weather Central's See dale. 


srid, with outline of North Pacific Ocean 
rectangular grid area used in study. See 
texte 
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Percent frequency of occurrence of fog in any form - 
1978). 


July (Guttman, 


Figure 4. 
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Percent frequency of occurrence of fog in any form - 
1978). 


August (Guttman, 


Figure 95. 
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OOOO GMT 3 August 1976, 


FNWC's Evaporative Heat Flux, EHF, 


North Pacific Ocean. 


Figure 7. 
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FNWC's Sensible/Evaporative Heat Flux, SEHF, 0000 GMT 3 August 1976, 


North Pacific Ocean. 


Figure 8. 
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0000 GMT 3 August 1976, 


ridional Wind Component, VY; 


ed) Me 
Ocean. 


FNWC's (deriv 


iene 1 


North Pacific 
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066° 


FNWC's Advective Fog Probability FTER, 0000 GMT 3 August 1976, 


North Pacific Ocean. 


Figure 10. 





‘o76T ysusny © LWNO 9000 ‘(ndvo 80L 0 


‘uvo0g otytoed UZION 
=Ap GSO) 4 aids Gage tase oe Ayrttaeqodd) I bvovod 7292 Ge lea 
OL6T qysn3ny woss poaAatsog uot yenby uotssei3say SdN wory APTTPAeqotd soq ‘tt ain3 TJ 
e.. ~ to” * =, 2 oe ; 


a 


SD, 


590 








‘or6T zsadsny € LN oo000 (CaHaS * 


Nic), Oat anes GOL’ S-GSEe 96 = AYTT TAS 
9161 psnsny WOTT pPoAT4og uot yenby uotssoid3y SdN wory AYTTTGUeqQoid 304 ‘St Site tat 
so 4 a hee 490 a a. “2D i ais 
7 : wae Op, 


repos 








Station 
Number 


32174 
32186 
32195 
32207 
62213 


* 
32215 


32217 
32559 
32594 
32618 
70454 


LAND STATIONS USED IN STUDY 
(U-5.uAIr Force. 


Latitude 
45.2 N 
46.2 N 
46.9 N 
43.3 N 
SIONS ae 
90.7 N 
00.0 N 
ae N 
Sino N 
54.9 N. 
o1.9 N 


APPENDIX A 


Longitude 


147. 
150. 
Sk 
153. 
156. 
156. 


Nnwo moO FPN NW DO OH © 
S= Aww ww HY Yt 


es 


Station not listed in reference. 
listed are as given with each synoptic report 
furnished by FNWC. 


no TZ) 


aK 


Surface 
Elevation (ft) 


38 
To 
26 
oI) 
42 
unknown 
ral 
88 
6 
ir 
14 


Lat/long values 


2 > 
Location given slightly removed from actual loca- 


tron OeooegeN Lata tide: 


accommodate study area. 


10) 


6620 ES vongitude to 





APPENDIX B 


OUTPUT PARAMETER DESCRIPTIONS 


Source Model 


Name 
Symbol Description Units 
A, Analysis Parameters (FNWC's Mass Structure Model) 
PS Sea Level Pressure (mb ) 
Analysis of observed sea level 
parameter. 
TAIR Surface Air Temperature Ce 


Analysis of observation-level 
air-temperature. 


EAIR Surface Vapor Pressure (mb ) 
Analysis of observation-level 
vapor pressure derived from the 
dew point. 


T925 925 mb Air Temperature (°C) 
Analysis of 925 mb air tempnera- 
ture. 

TSEA Sea Surface Temperature C205 


Once-daily analysis of observed 
sea-surface temperature. 


B. P,E. Parameters (FNWC's Primitive Equation Model) 


TX Surface Air Temperature Ge) 
Derived from surface air and 
potential temperatures, boundary 
layer depth, upper-level winds 
extrapolated to surface, air 
density, drag coefficient, gusti- 
ness factor, and empirical con- 
stants. 


EX Surface Vapor Pressure (mb ) 
Derived from model's mixing ratio. 


SOLARAD Solar Radiation (gm-cal/ 
Calculated absorotion of incoming 
Short-wave (solar) radiation 
(positive downward) 


om? hr) 
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fm . 


EHF 


SEHF 


SHF 


THE 


VVWW 


DDWW 


Spectral 
HW 

DPW 

PDW 

SPW 

SDW 


WCP 


Evaporative Heat Flux 

Derived using air density, drag 
coefficient, extrapolated winds, 
and mixing ratios. 


Sensible Plus Evaporative Heat Flux 
SEHF = SHF + EHF 


Sensible Heat Flux 

Recovered from SHF = SEHF — EHF. 
Originally derived by FNWC using 
drag coefficient, extrapolated 
winds, surface air temperature, TX, 
density, and constants. 


Total Heat Flux 

THF = SEHF - SOLARAD + LW, where 
LW is the heating due to long-wave 
(terrestrial) radiation. 


Marine Wind Model (FNWC) 


Marine Wind Speed 

Marine Wind Direction 

Both variables derived from a 
dynamic balancing of surface wind 
and sea-level pressure. 

Ocean Wave Model (S.O.W.M.)(CFNWC) 
Significant Wave Height 

Primary Wave Period 

Primary Wave Direction 

secondary Wave Period 


secondary Wave Direction 


Probability of White Caps 


Other Model Output Parameters ( FNWC ) 


SSTA 


Sea Surface Temperature Anomaly 
Calculated anomaly of sea-surface 
temperature from the mean of the 
day as interpolated from the 
monthly mean values. 
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(gm-—cal / 


ene hr ) 


(gm-cal / 


om? hr) 


(gm-cal/ 


eine hr) 


(gm—-cal/ 


em Hr) 


(knots) 


(degrees / 
10) 


(feet) 
(sec) 
(degrees / 
10) 
(sec) 
(degrees / 


iO) ) 
(percent ) 


Ge 





Derived Parameters 


u Zonal Wind Component (m/sec) 
u = -VVWW sin(DDWW = 10) 


V Meridional Wind Component (m/sec) 
v = -VVWW cos(DDWW - 10) 


CAPU I Directional Wind Component (m/sec) 
CAPU = -u * sin(LNGA) - v + cos(LNGA), 
(Hadtaner. 19/1). 


CAPV J Directional Wind Component (m/sec) 
CAPV = u « cos(LNGA) - v « sin(LNGA), 
(Haltiner, 1971), where LNGA = 
-10 - (I,J point longitude). 


THETAX Potential Temperature X C8) 
Derived using PS, TX. 


THETAR Potential Temperature R Gie 
Derived using PS, TAIR. 


STABX Stability xX (°K/mb) 
Derived using ({THETAX-(THETA of T925) ] 
/(PS-925). Value greater than zero 
indicates absolute instability. 


STABR Stability R (°K/mb) 
Derived using (THETAR-(THETA of T925) ] 
/(PS-925). Same value effect as STABX. 


ASTDX Air-Sea Temperature Difference X eee) 
ASTDX = TX — TSEA 

ASTDR Air-Sea Temperature Difference R Ge 
ASTDR = TAIR —- TSEA, 

ADTSEA Advection of TSEA CAC 
Formulae and notes below. Hour ) 

ADTX Advection of TX ney) 
Formulae and notes below. Gur ) 

ADTAIR Advection of TAIR. Go) 
Formulae and notes below. HOW ) 

AASTDX  Advection of ASTDX Cae 
Formulae and notes below. Hour ) 


G3 





AASTDR Advection of ASTDR CH) 
Formulae and notes below. Hour ) 
G. Climatological Parameters 
CLIMO National Climatic Center Fog (percent/ 
Frequency Climatology 100) 
PERS Persistence of Fog Likelihood (percent / 
100) 
H. Interactive Parameters 
CLISEHF CLIMO «+ SEHF (% gm-cal/ 
100 om“hr) 
Ghiy CLIMO + v (m %/ 
100 sec) 


VSEHF =v -: SEHF £100 em-cal/ 
So ee = 


em sec hr) 


SEHFSQ SEHF - SEHF (gm -cal / 
em” br) 
CLIEHF CLIMO - EHF (% gm-cal/ 
100 aa eb 
CLIADT CLIMO - ADTAIR (% -C/ 
100 hr) 
EHFSQ EHF - EHF (gm -cal / 
yal bisa) 
EHFADT EHF - ADTAIR (gm-cal °C/ 
em” hr) 


Advection Formulae and Conditions: 
For the advection of a quantity (R) the formula, ADQ = 


-V-V(Q), was used in the following finite difference form: 


_ _ RMAP 
cee nee sree te ee 854) ~ 2 ySq) 7) 





where RMAP 


tl 


(1l+sin(60))/(C1l+sin( latitude) ) 


sand DM = [(2)*(6.37*10°)+*(1+sin(60))]/31.205 


(31.205 = grid mesh lengths, pole to equator, on FNWC's 
ie ieorid ). 


In the temperature advection calculation for point C, using 
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the five grid points illustrated below, one or two of the 
points, namely, A, B, D, or E, may be outside the study area. 
In the bogusing method suggested by Mr. Leo Clarke, FNWC, 
When a non-center point (e.g., point A) was judged to prob- 
ably produce a land and/or dissimilar sea area influence on 
the resulting advection, the center point value (point C) 
was substituted for it. This "bogussing" is necessary to 
maintain a "purely'' marine characteristic in the resultant 
parameter value. 

The set of study area boundary grid points used for 
bogusing (and some double bogusing) are depicted in Fig. l. 
The upper case letters denote the points whose values were 
used for bogusing. The lower case letters mark the posi- 
tions of the adjacent points being bogused. The only bound- 
ary point close enough to land to give concern is the one 


near station number 32594 (south tip, Kamchatka Peninsula). 
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APPENDIX C 


Abridged version of internationally used weather code 
figures and definitions for reporting present and past 
weather and low clouds in the surface synoptic report 
(U. S. Departments of Commerce, Defense, and Transporta- 
tion, 1969). 


Present Weather 


Code 

Value Definition 

00-03 Characteristic change of the state of the 
sky (cloud) during the past hour. 

04-09 Haze, dust, sand, or smoke. 

10 Deep light fog. 

11-12 Shallow heavy fog. 

f3—) 7 Lightning, thunder, or precipitation within 
Sachin OL reaching fhe sround. 

18-19 Squall(s), funnel cloud(s) during the past 
hour. 

20 Drizzle during the past hour. 

21-23 Rain, snow, or rain and snow during the past 
heune 

24 Peeezine Grigzte GQuring the past hour. 

29-27 Shower(s) during the preceding hour. 

28 Fog during the past hour. 

29 Thunderstorm during the past hour. 

30-39 Duststorm, sandstorm, drifting or blowing snow. 

40 Fog at distance, but not at station, during 
the past hour (visibility less than 1 km). 

41-49 Deep heavy fog at the time of observation 
(visibility less than 1 km). 

20-59 BbElZzZzre.  Oruarl az ewan radn . 

60-63 Slight to moderate rain. 

64-65 Heavy rain. 

66 Slight freeZing rain. 

67 Moderate or heavy freezing rain. 

68 Slee cat or wart 27 ee land- snow. 

69 Moderate or heavy rain or drizzle and snow. 


76 





70-~—79 Solid precipitation not in showers. 


80-99 Showery precipitation or precipitation with 
current or recent thunderstorms. 


Past Weather 


0 Cloud covering 4 or less of sky throughout the 
period. 

IL Cloud covering more than $ of sky during part of 
the period. 

2 Cloud covering more than 4 of sky throughout period. 

3 Sandstorm, or duststorm, or blowing snow. 

4 Heavy fog, thick haze, or smoke. 

3) Drizzle 

6 Rain 

7 Snow, rain and snow. 

8 Shower(s). 

9 Thunderstorm, with or without precipitation. 


Low Cloud Type 


O No low clouds. 
a Ragged cumulus of fair weather. 
2 Generally towering cumulus. 
3 Cumulonimbus without cirriform or anvil tops. 
4 Stratocumulus formed by cumulus spreading out. 
5 Stratocumulus not formed by cumulus spreading. 
6 - Stratus or fractostratus. 
a Fractostratus of bad weather. 
8 Cumulus and stratocumulus, with bases at different 
levels. 
g Cunulenambus with clirriform top. 
/ Low cloud obscured. 
fiei bility 


20 Less than 50 m 

oa. O-199 m 

92 200-499 m 

93 200 m —- 0.99 km 

94 w=. Oo km 

2s 2—- 3.99 km 

96 4 —- 9.99 km 
97-99 equal to or greater than 10 km 
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APPENDIX D 


FOGCAT I Categorization Scheme 
(Van Orman and Renard, 1977) 


1. Groupings and symbols used in FOGCAT I categorization 


scheme. 
Present Weather (ww) Past Weather (W) Low Cloud (CL) 
Symbol Associated Symbol Associated Symbol Associated 
ww Codes W Codes CL Codes 
41G 41-49 4 4 6 6 
10G 10,28,40 4,95 4,9 Out Sirti 
1G ieee. 20, 24 5G OMS B i 
0G 5D0-59 2G On 2 7 any Cl 10 
listed above 
60G 60-63,66,68 a: any W not 
listed above 
“0 any ww not 


listed above 
G = Groun B = Low clouds obscured 
2. 


Scheme for categorizing observations according to like- 
lihood of fog (FOGCAT I). 


Major Sub- Present Past Low Cloud 
Category Category Weather Weather Type 
(ww ) (W) (CL) 
Strong al 41G 4,5 6 
Foggers a2 ss . B 
= § bl ‘A 2G 6 
b2 tt tt B 
b3 10G 4,95 6 
D4 tt tf B 
el 41G * 6 
2 A 4,5 Sy 5 we 
c3 tt Tf * 
c4 10G 2G 6 
Foggers dl 41G “s B 
oe d2 7 2G a, Th 
d3 o0G 4 6 
d4 1? tt B 
d5 10G 3 6 
d6 is 2G B 
d7 a 4 ene 


~“e 
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el 41G 
tt 


e2 
e3 290G 
e4 10G 
e5 tt 
e6 L1G 
a7 tt 
Siig 41G 
f2 20G 
2 LOG 
£4 11G 
To tt 
gl LOG 
g2 tt 
g3 11G 
o4 tt 
g5 t 
hl tt 
h2 Tf 
geil 19G 
12 11G 
svral tt 
Past/Weak kl 90G 
Foggers k2 60G 
= P v1 S 
ce = 
m1 20G 
m2 = 
nl = 
el 60G 
o2 * 
pl 60G 
p2 es 
p3 is 
Maybe ql 20G 
Foggers q2 i 
= M q3 tt 
q4 tT 
qo 11G 
q6 60G 
Giz tt 
ik 290G 
r?2 tt 
r3 60G 
foce tt 
r5 tt 
r6 tt 


ION 


“ee 


~] 


“ee 


~] 


Ne 


~] 


“oe 


~] 


“ee 


~] 


“e 


¥eOKOWD KOH EOI MDNHKK*#WH*EWDN KH 


“eo “eo 


EHD *F*UNUAHWANAHWHAD 
IN 


PPP HOR * 


IN 


“ee “e 


w 


WOW *WAM * xO 
ION 


“ee 





Non- el: 60G OG * 

Foggers u2 * = 

= N vil * OG B 
wl * B 
x1 . a5 2). Tf 
yl 2 OG Seer 
y2 * 5G * 
Z1 > * *K 


3. Major fog categories and sub-category fog groups with 
associated Fog Probabilities as a function of 
Vis LD mie iy. 


Visibility Code Values 


96-99 94-95 90-93 Fog | 
(Good) Giiavine) (Poor) Probaba bit y 


Major 
Category Subcategory Groups 


not 
used 


+ulyols on Ben br oa © ale Oo py 





nine cele aly o bo om Bm ban rm no alo oe 
on 
on 


NS wade 
Hs 
ca 
OROWANNAKRFPOAOUOAWWANAKPUODOAROWWNARAO 


AS. 
pe wade [ealrob opaerbrraned oe 
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APPENDIX E 


FOGCAT II Categorization Scheme (Quinn, 1978) 
Fog 
Present Weather Past Weather Valsacioaiad & y Probability 

ww ) (CW VIS i) 
0-9 - _ Q 
10-12 - ~ 100 
13-19 - — 0 
20 0-3 ,6-9 - 0 
20 4,5 90-95 S15) 
20 4,9 96-99 0 
21-24 0-3 ,6-9 _ 15 
21-24 4,9 90-95 60 
21-24 4,5 96-99 iS 
25-27 0-3,6-9 — 0 
25-27 4,5 90-95 10 
25-27 4,95 96-99 0 
28 ~ - 100 
29-39 — - 0 
40-49 — = 100 
90-59 0-3,6-9 = Sa 
50-59 4,5 90-95 85 
50-59 4,5 96-99 35 
60-69 0-3,6-9 - 15 
60-69 4,5 90-95 60 
60-69 4,5 96-99 sis) 
710-79 _ - O 
80-89 0-3 ,6-9 ~ 0 
80-39 4,5 90-95 10 
80-39 4,9 96-99 0) 
90-99 = ~ 0 


Dash indicates that the particular category was not con- 
Sidered in the assignment of a fog probability. 
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APPENDIX F 


Verification Score Formulae (Quinn, 1978) 


EVENT ESTIMATED 


YES NO 
YES A g Total (T) = A+B+C+D 
EVENT 
OBSERVED No. of Correct 
NO B D Forecast FC = A+D 
Heidke - FC — EX , 7 2BC 
Skill Score ~T— EX Been mee ee 


where EX = (AFB) CAtC) PSB Or) ? 


No. of expected 
correct forecasts due to chance. 

ee) ee 

T-D A+B+C 


Both scores indicate more skill with larger 
positive values. 


Threat Score = hance: 005 15 5 2 


The Probability Score (PS) is from that given by 
Panofsky and Brier (1958) and may be written as (Renard, 
EO 7a) 


gl 


nN 
O 

2 2 2 
shoal 


‘oe 


P-—Score 


Ci. 
I 
j! 


Range: O < PS < 2 
where 
N = Total number of cases 
n = Total number of non-events 
Pp. = Associated probability value for the non-event 
n, = Total number of events 
p. = Associated probability value for the event 


The closer to zero, the greater the skill. 
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